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For Rayleigh-Benard convection of a fluid with Prandtl number er sa 1, we report experimental and 
theoretical results on a pattern selection mechanism for cell-filling, giant, rotating spirals. We show 
that the pattern selection in a certain limit can be explained quantitatively by a phase-diffusion 
mechanism. This mechanism for pattern selection is very different from that for spirals in excitable 
media. 
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Spiral patterns are found in many pattern-forming sys¬ 
tems jjj. Famous examples include cardiac arrhythmias 
of the heart Q, the aggregation of starving slime mold 
amoebae |jj, and the Belousov-Zhabotinsky chemical re¬ 
action Many of these systems can be classified as 
excitable media in which the core of the spiral, like a 
pacemaker, selects the temporal and spatial evolution of 
the outward travelling spiral waves Q] . In this Letter we 
present a detailed experimental study of a driven, dissi¬ 
pative system in which the formation of spirals can be 
attributed to a qualitatively different mechanism acting 
far away from the spiral’s core. In particular, we show 
for Rayleigh-Benard convection (RBC) of a small Prandtl 
number fluid that the rotation of giant, multi-armed spi¬ 
rals can be captured using concepts based on nonlinear 
phase equations IB- These concepts should be uni¬ 
versal and preliminary evidence indicates that similar 
reasoning may also apply to the spiral pattern formation 
in vibrating granular layers 0 and in gas discharges 

0 - „ 

As stated earlier Jl2[ , the rigid rotation of a giant, fi¬ 
nite spiral of radius rd necessitates that the spiral waves 
which propagate from the spiral’s core are annihilated at 
r = ra by a circular motion of the outer defect. Thus, the 
pattern simply consists of stationary, concentric rolls for 
r > rd- This balancing mechanism has been placed on a 
more precise theoretical framework by Cross and Tu [|7l|| 
(CTC). They argue that the rotation of a spiral requires 
the reconciliation of two competing selection principles 
acting far away from the spiral’s core: (1) wavelength se¬ 
lection by climb of the outer defect and (2) the emission 
of radially travelling waves due to target selection. These 
arguments have been successfully tested within general¬ 
ized Swift-Hohenberg (SH) models for axisymmetric spi¬ 
rals |t] ^.[bi|l and have also been extended to multi-armed 
spirals |)|. However, the generalized SH-model used in 
these studies contains ad hoc parameters, and its stabil¬ 
ity regime deviates substantially from that for RBC [[wj. 
As Cross |U has pointed out in his concluding remarks, 


serious conceptual uncertainties remain as well. 

In this Letter, we present the first quantitative exper¬ 
imental and theoretical analysis of giant rotating spi¬ 
rals in RBC. For axisymmetric, multi-armed spirals in 
large aspect-ratio cells we find good quantitative agree¬ 
ment between experimental measurements and theoret¬ 
ical predictions. However, for the frequent case of non- 
axisymmetric spirals our measurements are in conflict 
with the proposed target selection mechanism and, as a 
consequence, CTC’s “Invasive Chaos” idea as a tentative 
explanation of Spiral Defect Chaos (SDC) |h| requires 
refinement. 



FIG. 1. Shadowgraph image of the six convection cells for 
t = 0.98. Dark corresponds to warm up-flow, while light 
corresponds to cold down-flow. The cells used in the analysis 
are numbered from 1 to 4. Cell 3 contains a PanAm pattern 
and cell 2 contains Spiral Defect Chaos. 

Rayleigh-Benard convection occurs in a horizontal fluid 
layer of height d heated from below when the temper¬ 
ature difference AT exceeds a critical value AT C . For 
e = (AT/AT C — 1) > 0, a pattern of convection rolls with 
wavenumber k « n/d develops [[HJ. Recent large aspect- 
ratio experiments using pressurized gases revealed inter¬ 
esting nonrelaxational pattern evolution. For moderate e 
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Bodenschatz et al. |12| found giant rotating spirals simi¬ 
lar to those described here, while Morris et al. jL5| found 
the spatiotemporal chaotic state of spiral defect chaos 
(SDC). 

As shown in Fig. [I], our experiment consisted of six 
circular convection cells filled with pressurized CO 2 gas. 
The convective pattern was visualized with the shadow¬ 
graph technique [^6|. The experimental setup including 
parameters is described in detail in Ref. 1171. In most 
experimental runs reported here we increased e quasi- 
statically starting from below onset of convection (e = 0). 
Upon crossing onset, a small amplitude axisymmetric 
target patterns developed in all cells due to weak static 
sidewall forcing. (We achieved this by a step in the cells’ 
sidewalls JT^]). Above e « 0.4 the initially axisymmetric 
target moved off-center, compressing the pattern on one 
side while dilating it on the other. The appearance of 
this instability is consistent with theoretical predictions 
||. Upon further increases in e the target moved far¬ 
ther off-center until the wavenumber in the compressed 
region increased beyond the skewed-varicose (SV) insta¬ 
bility Q and a defect pair would nucleate to decrease the 
wavenumber. One of the defects then moved to the cen¬ 
ter while the other glided radially outward before coming 
to rest at a distance from the geometrical center. The 
pattern then relaxed to an on-center, one-armed rotating 
spiral of radius r^. (While this instability was observed 
for e « 0.5 in the largest cell, it was slightly postponed 
in the smaller cells). Above e « 0.55 we observed a be¬ 
havior reminiscent of the target instability in which the 
spiral’s core would move off-center. For e < 0.64 we 
observed stable, rotating, off-center spirals (see Fig. 3 
below) which with each increase in e would move fur¬ 
ther off-center. Eventually, the wavenumber in the com¬ 
pressed region increased beyond the SV-instability and 
defect pairs nucleated. The pattern then developed into 
a three- or four-armed spiral, into the so-called PanAm 
pattern, or into Spiral Defect Chaos. Examples of the 
latter two are shown in Fig. |l| in cell 3 and cell 2, re¬ 
spectively. We note that stability properties of spirals 
had not been addressed experimentally nor theoretically 
prior to our investigations. 

For a few additional runs, we jumped the control pa¬ 
rameter from below the onset of convection (e < 0) to 
above (e > 0). As the jump was increased, we observed 
targets, one-armed spirals, multi-armed spirals, PanAm 
patterns, and SDC. Interestingly, we observed two-armed 
spirals only when we employed this procedure. Other¬ 
wise, the general trend with e observed using the two 
methods was similar. We note that when e was quasi- 
statically decreased starting from SDC or PanAm pat¬ 
terns, we observed PanAm patterns, not targets, even 
close to onset. This shows that the static sidewall forc¬ 
ing was very weak. 

Let us initially focus on the one-armed spiral. Fig. |2j 
shows the average velocities Vd of the outer defects for 


an experimental run with one-armed spirals in cells 1-4. 
(The velocities are normalized in terms of cL/tt, where tt 
is the vertical thermal diffusion timescale.) Note that the 
spirals had a variety of different sizes r^. We measured 
Vd by tracking the path of the outer defect with a Fourier 
demodulation technique [|I8|,^9|. For on-center spirals the 
averaged defect velocities of spirals in cells 1-3 obeyed 
the same linear relationship Vd = 0.64(e — 0.09) . For 
values of e with off-center spirals, Vd changed abruptly. 
(The spiral in cell 4 showed deviations from this behavior, 
possibly due to its small size.) We also measured the 
rotation frequency ui of the spirals and found that the 
data are well-described by u> = Vd/rd■ This dependence 
on spiral size differs markedly from the size-independence 
of uj observed for spirals in excitable media 0. 



FIG. 2. Average velocity Vd of the outer defect vs. t for 
one-armed spirals. The size of the on-center spirals is given 
in Ref. |^|. The dashed line represents a linear fit to the 
on-center data. 

The first part of CTC’s argument is that the behavior 
of the outer defect of a spiral can be considered as a dis¬ 
location climbing in a (slightly curved) roll pattern with 
a radial wavenumber q(r) in analogy to the simpler case 
of a dislocation climbing with the velocity Vd in a straight 
roll pattern with background wavenumber q. This latter 
case was considered theoretically in the framework of the 
Swift-Hohenberg model ||. Based on the phase diffusion 
equation which captures the behavior far from the defect, 
the defect velocity Vd was found to obey the relation |^2) : 

Vd = P{e)(q - Qd(e)), (1) 

where qd is the zero-velocity wavenumber. 

In order to compare with this theoretical prediction 
one first needs a sensible definition of a background 
wavenumber q(r). Crucial to this was the observation 
that the time-average of an m-armed spiral (when per¬ 
formed over a duration equal to a multiple of the ro¬ 
tation period) yielded a target pattern. This is exem¬ 
plified in Fig. 0(B) for the one-armed off-center spiral 
shown in Fig. 0(A). As shown below, choosing the radial 
wavenumber q(i ■) of the target pattern as the background 
wavenumber q in Eq. ([!]) appeared to work quite well. 

The occurrence of a target pattern after averaging over 
a rotating spiral can be rationalized by the following ap¬ 
proximate calculation. A one-armed, finite spiral can be 
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described satisfactorily by a modified Archimedean spi¬ 
ral, i.e., f(r,t) = A(r) cos (qr — <j> + 4>d ~ wi), where q 
is the wavenumber of the underlying target without de¬ 
fects, A(r) is the amplitude, and r is the radial distance 
from the spiral’s center. The phases tf> = arctan(y/a;) and 
4>d = arctan((y — yd)/(x — Xd) are polar angles centered 
about the spiral’s core at r = (0, 0) and the outer defect 
position fd(t) = ( Xd,yd ), respectively. It is not difficult 
to show that the time-average of such a rigidly rotating 
spiral gives J(r, r^) = 27rA(r)Ji (r/rd) cos(gr), where J\ 
is the first Bessel function of the first kind. We note 
that the corresponding average of an Archimedean spiral 
{<j>d = 0) would vanish. 


FIG. 3. (A) Off-center one-armed spiral and (B) the av¬ 
erage over one rotation period overlaid with the defect posi¬ 
tions at constant time intervals At = 8.7 tt for e = 0.62 and 
a = 1.38. 

Using the above definition of the background 
wavenumber q(r), we tested Eq. ( 0 ) using experimen¬ 
tal velocities obtained from off-center, one-armed spirals 
over a wide range of values of e. The slower defect mo¬ 
tion in the dilated regions and the faster motion in the 
compressed regions seen in Fig. 3B is consistent with 
Eq. ([j]). For a quantitative test we measured q(r) at the 
defect positions by fitting small regions to concentric roll 
patterns 0. (The method of Ref. 0] produced simi¬ 
lar values.) We then plotted Vd versus q for each value 
of e and found that the data were well-described by a 
linear relationship, allowing us to determine the parame¬ 
ters /3(e) and qd(e) using a least-squares fitting procedure 
fl4| . These data are shown in Fig. [|. Using simulations of 
the three-dimensional Boussinesq equations , we also 
measured defect velocities Vd as a function of background 
wavenumber q for defects moving in straight roll patterns. 
Again, we found a linear relationship between Vd and q 
and determined /3(e) and qd(e). As shown in Fig. |], the 
numerically determined values for /3(e) and qd(e) are in 
excellent agreement with the experiment. Using the sim¬ 
ulations, we also analyzed off-center, one-armed spirals 
pb[ and found excellent agreement between experiment 
and theory Q . 

The second part of CTC’s argument — selection by cir¬ 
cular traveling waves — relies on the fact that away from 
the core of the spiral the wavefronts deviate only slightly 
from circular and are thus well-approximated as targets. 
It has been shown that targets prefer a specific wavenum¬ 


ber q t (e) and that a target with a wavenumber dif¬ 
fering from qt(e) will attempt to adjust its wavenumber 
by emitting circular waves of frequency w t (e). Using the 
nonlinear Cross-Newell phase-diffusion equation ||, one 
finds 

u t = a(e)(q t (e) - q)/r, (2) 

where a = 2D\\(q t ), D\\(q) is the parallel diffusion con¬ 
stant, and r is the distance from the center of the target. 
The numerical value of the parameter a can be calculated 
from the growth rate cr(q t , K) of a longitudinal modula¬ 
tion with wavenumber K of a pattern with wavenumber 
qt as a = (~2a(q t , K)/K 2 )\ k ^ 0 [|]. 



FIG. 4. (A) (3 and (B) qd vs. e for single-armed, off-center 

spirals, compared with numerical simulations of defects in a 
straight roll pattern. 

Using geometrical arguments, CTC showed that a one- 
armed spiral requires: 

Vd = v t r d , (3) 

where, again, Vd is the velocity of the defect at radius 
rd- We used Eq. 1) and Eq. (|2j) and the generalization 
to multi-armed spirals jfij to determine a. From the in¬ 
vestigations above and a similar analysis for m- armed, 
on-center spirals, we measured the average defect veloc¬ 
ity Vd- For each cell we determined qt by extrapolating 
the azimuthal average of qt (e) at the defect positions from 
on-center targets at lower e to the larger e where we ob¬ 
served the m-armed, on-center spirals. We found that 
(Qt - q c ) ~ e The data for q t was in good agreement 
with the numerical predictions by Buell and Catton [^8) . 
The remaining unknown background wavenumber q was 
measured from the underlying target, again azimuthally 
averaged over the defect positions. With this information 
we determined a for one-, two-, three-, and four-armed 
on-center spirals. The data are summarized in Fig. |. It 
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shows that for the larger cells and the multi-armed spirals 
a agrees well with the theoretical value numerically de¬ 
termined from the growth rate of the Eckhaus instability 
using a Galerkin code (dashed line). The discrepancies 
become significant for one-armed spirals in the smaller 
cells 3 and 4. 



FIG. 5. a vs. t for single- and multi-armed spirals. 
The shapes of the symbols correspond to the four different 
cells. The level of gray in the symbols distinguishes between 
the number of arms in the spirals. One-armed spirals are 
clear, two-armed are gray, three-armed are dark gray, and 
four-armed are black. The numerically determined value of a 
is plotted as a dashed line. 

For the frequently occurring off-center spirals we found 
that the rotating spirals had regions where the local 
wavenumber was larger than that of the stationary tar¬ 
get. This observation indicates that a “local” interpreta¬ 
tion of Eq. (2) (negative a; t ) is not reasonable. Rotating 
spirals are important structure elements of SDC ||l5| . ^3| . 
Mostly they are one-armed and often strongly off-center. 
In this case the target selection becomes unreliable and 
as a consequence CTC’s “Invasive Chaos” idea as an ex¬ 
planation of SDC has to be reconsidered. 

In summary, our analysis presents the first quantita¬ 
tive test of Eq. ( 0 ) 1 and Eq. ( 0 ) 01 in a real physical 
system. Giant multi-armed spirals are exceptionally well- 
suited for the study of dislocation dynamics in RBC; one 
can follow the defect trajectories for very long times and, 
in addition, the background wavenumber is naturally de¬ 
fined by the time average of the spiral. 
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